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The optimal use of the continuous review model requires 
that an inventory system be examined after the receipt of 
every demand. Sometimes a delay can be encountered in the 
timing of these reviews due to several uncontrollable 
factors. As the length of these delays increases, a point 
1s reached where it 1s better to switch the inventory system 
to a periodic review model. This thesis develops a methodology 
by which this point can be found for varying cost factors and 
demand levels. Using simulation, an example series of curves 
is presented that demonstrates the optimal point to switch 
inventory models for selected lead times. If delays are 
expected in tne time between reviews, using the methodology 
offered in this paper will provide the manager an 
informational criterion for deciding what inventory model to 
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A. BACKGROUND 

Traditional inventory systems revolve around two basic 
models. These models are the continuous review model, often 
called the Q-model, and the fixed period review model, often 
referred to as the P-model. Most large inventory systems 
use one of the two aforementioned models, in full or in part, 
to control inventories wnile minimizing the variable order 
and holding costs while also maintaining some minimum level 
of service. 

The continuous review model is so named because inventory 
levels are reviewed every time an item is demanded. This 
model allows a manager to be more aware of a potential 
out-of-stock condition, but it can require a costly review 
system. Ina fixed period review model inventory levels 
are checked only after a set period of time has elapsed. 

The total yearly costs of inventory review is less than what 
would be required for the Q-system, but large variances in 
demand could cause additional stock-outs if the time between 
reviews is extensive. 

With the advent of the computer more and more activities 
use the continuous review model [Ref. 1]. They feel that the 
additional responsiveness in the Q-model compensates them 


(in reduced stock-outs) for the more expensive review system. 





Tnis trade-off between the implicit costs of a stock-out and 
me actual costs OE the review can be easily represented, 
and quantified, with optimal cost formulas theoretically 
determined oe each model. However, few systems are operated 
in a totally optimal manner. Most organizations find them- 
selves constrained by uncontrollable variables and conflicting 
objectives that require thelr inventory systems to be operated 
in a less than optimal fashion. Each one of these additional 
constraints exacts a cost that must be calculated in order to 
determine the true total costs of operating an inventory 
system. The purpose of this thesis is to examine one of 
tnese non-optimal inventory systems and show, through simula- 
tion, what additional costs result when the system is not 
operated in a totally optimal manner. 

The United States Navy uses a continuous review system 
at their Inventory Control Points (ICP). Theoretically, the 
inventory level of each stocked item is checked upon every 
issue. If an item is below a pre-determined reorder point a 
requisition is generated to the appropriate procurement 
activity which, at some later date, will bring the inventory 
back up above the reorder point [Ref. 2]. Many times, 
however, an individual item's inventory is not examined when 
there has been demand activity. Several factors can cause 
delays in these examinations. Some of these factors are 
computer downtime, delays in funding, or inadequate personnel 


resources to do the reviews. When delays occur between 








receipt of a demand and a check of the relevant inventory 
level the continuous review system is not being operated 


optimally and a cost is being incurred by the Navy. 


Bae URPOSE 

This thesis examines the costs of running a continuous 
review model non-optimally by suing simulation to create an 
inventory system similar to the Navy's. The first simulation 
is configured to show both the theoretical optimal (lowest 
cost) inventory system, and various ee eclons where there 
are fixed period of review imposed on the continuous review 
model. Since it is possible that the period between reviews 
could extend long enough in the Q-model that the fixed period 
review model would become more applicable to a given situa- 
tion, a second simulation model was developed. That model 
imitates wnat a Navy inventory system would be like if a 
Pemodel was employed at the Inventory Control Points for 
specific items. The use of these two simulations gives us 
the ability to compare and contrast the Q and P models in 
order to determine which model is best for any given item 
and period of review. In the next chapter, the specific 
make-up of each one of these models is outlined and an 
analysis 1s completed on the mathematical differences 
between the Q and P system as well as the inherent advantages 
and disadvantages of each model. Chapter III describes what 


assumptions were used in the development of the simulations. 





Appendices A and B present the programs used to make the 
Simulation runs. 

The last three chapters provide an analysis of the types 
of information available from the simulations. Chapter IV 
details a methodology for constructing a curve comparison 
between the two models using the simulation programs. The 
following chapter then takes select data from these curves 
and presents an analysis of trends, optimal model costs, and 
costs to do a review. The last chapter.contains a summary 


of the thesis along with conclusions and recommendations. 
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Pa ONT RODUCTION 
In general, inventory theory concerns itself with two 

models--the continuous review model (Q-model) and the-fixed 
period review model (P-model). These two models can be 
considered under either assumptions of deterministic or 
stochastic demand. Since the purpose of this thesis is to 
Simulate, as close as possible, an actual inventory system, 
the version having stochastic demand will be used. The 
associated mathematical theory and imputed optimal cost 
calculations for each model are outlined in the sections 


below. 


Eee ONTINUOUS REVIEW MODEL 

The basic element common to continuous review models is 
that the state of the supply system is known at any instant 
of time and, therefore, decisions may be made at the instant 
that an inventory position (on-hand inventory plus on-order 
minus backorders) falls to some reorder level. In the 
Q-model an optimal quantity Q, based on the Wilson Economic 
Order Quantity (E0Q), is assumed to be ordered whenever the 
inventory position falls to a reorder level (r). Thus the 
inventory position is raised to Q + r immediately following 


imme order [Ref. 3]. 


ieee 





When working with the stochastic Q-model, the inventory 
can be divided into working stock and safety stock. Working 
stock is what is expected to be used during a given time 
period. The average working stock is one-half the optimal 
order quantity. Safety stock does not depend on the optimal 
order quantity. It is based on the variability of the 
stochastic demand distribution and is used to protect against 
higher than expected demand levels. 

Each unit of increase in safety stock provides a demin- 
ishing benefit. The first unit of inventory above the 
expected demand provides the largest marginal increment of 
protection against running out of stock. As additional units 
of safety stock are added this marginal protection gradually 
decreases until the cost of storing additional units plus 
the expected stock-out cost 1S at a minimum. This level is 
defined to be the optimal safety stock level and movement 
away from this result produces an increase in variable 
costs [Ref. 4]. 

In continuous review models the reorder point is composed 
of a mean leadtime demand (the time between placing an order 
and wnen the order arrives) plus the safety stock. Through- 
out this thesis the mean leadtime demand will be specified 
and a Poisson distribution will be assumed for the random 
demand [LRef. 5]. The Poisson probability distribution was 
selected over a normal probability distribution because; 


(1) it more accurately represents items with low demand 
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(less than 20/year), (2) at higher demand levels the Poisson 
approximates the normal even though some round-off problems 
can occur, and (3) the Poisson avoids the possibility of a 
negative demand. 

The actual determination of the optimal order quantity 
(Q) and reorder point (r) requries an iterative approach 
because of the discrete demand distribution. The optimal 


values are the largest Q such that [Ref. 6]: 


2.1 [A + mn(r)] 


OCQ-1) < ic 


and the largest r such that: 


mcr) > STOP ay 
where: 
T = stock-out cost per unit short 
r = reorder point 
A = order cost per order 
C = purchase price of one unit 
IC = holding cost per unit 
Q = opermal order quantity 
d = expected annual demand 


H(r)= the probability of a stock-out during a 
procurement lead time when r is the reorder 
DOr e 


n(r)= expected number of stock-outs per order cycle 
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Mtessteps Of the approach are: 





1) Find largest Q such that: Q(Q-1) < ate 
ee trimeric Oldetermine RISK = ———— 
: CLG 7h 


Using this RISK value determine the largest value of r 


auen that: 


R 


CCDF = p(x) > RISK 


uo 


gy 


where p(x) is the probability of a demand of x during 
procurement lead time. 

A Poisson distribution table is needed because of the 
iterative process. When this table is generated (using the 
applicable mean lead time demand) it 1s necessary to 
establish the complementary cumulative distribution (CCDF) 
Of the distribution. Appendix III is a program that will 
compute the Poisson CCDF using different mean leadtimes. 

3) Using the above r, and the CCDF Poisson ela ieycicalleybles eg 


mame find: [Ref. 5] 


mee) =<9(M = CCDF(r=1)) = (r * CCDF(r)) 
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fas, 


(eewtthn this value of n(r) we determine the largest Q 


such that: 
2x {A + amn(r)] 
Q(Q-1) = oie 


5) Returning to step 2 with this value of Q, the steps 
are repeated until such time that the r does not Enotes with 
the next eee eton. The resulting Q and r are optimal and 
the safety stock can be computed by subtracting the expected 
mean lead time demand (M) from r. 


6) Compute the expected annual holding cost: 
holding cost = IC [3 + (r - M) + n(r)] 

7) Compute the expected annual order cost: 
order cost = on 

8) Determine the expected annual stock-out costs: 
stock-out cost = - 7 n(r) 


9) Add up the order, holding, and stock-out costs from 


oeep> G6, /, and 8 to get total variable costs. 


tTyvc = A484 Icle + Cee ee eee or ce) 


Q Q 
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The above total optimal cost information will provide the 
baseline costs by which comparisons will be made when 
changes are introduced into the simulation. 

Since the purpose of this thesis is to compare the 
P-model with the Q-model (constrained to operate as the 
P-model), the costs of doing the review must be considered. 
It 1s the product of the expected annual demand i and the 


eo eeere) Of an individual review. 


Review Cost = ARC 


However, review costs need not be considered as a direct 
part of the total variable cost equation in the analysis 
in Chapter IV . Their effect will be analyzed in Chapter \V. 
The Q-model iS more expensive to administer than the 
P-model. A perpetual and continuous inventory record must 
be maintained which often requires computer support for a 
large number of product items. However, this perpetual 
inventory system is the most demand responsive and this 
characteristic becomes even more important as demand 
variation increases [Ref. 7]. The Q-model also has a lower 
average inventory level because less safety stock is 
required than the P-model. An additional advantage of the 
continuous review model is that it keeps management in close 
touch with changing demand patterns which aide in production 


planning as well as in inventory management. 
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fee LXED PERIOD REVIEW MODEL 

Time based inventory models are most appropriate when an 
organization is faced with fairly uniform independent demand 
and wants to maintain an inventory most of the time. When 
demand is continuous time and quantity can be considered 
interchangeable [Ref. 4]. 

For the fixed period review model we need to calculate 
two variables--the time between reviews (T) and the maximum 
iiemeory level (FE). Only after the time T has passed will 
inventory levels be checked and an order placed. The size 
of the order will be the difference between the maximum 
inventory level for each item and the inventory position at 
the time of the review. Unlike the Q-model where the order 
quantity is constant, the order quantity for the P-model will 
vary depending on fluctuations in demand between orders. 

The maximum inventory (E) will consist of the safety 
stock plus the average demand during an order interval and 
lead time. In the fixed order interval system safety stock 
is needed for both the order lead time (M) and order interval 
@meetref. 4]. 

The Poisson distribution is also used in the fixed period 
review model to simulate random demand. However, the mean 
of the distribution represents demand both during the order 


eyecle and lead time. 


ey 


rat 





The steps for finding the optimal (T) and (E) also require 


mete rative solution similar to that of the continuous 


review model. The following is the solution procedure for 


the P-model 


Bret. out 


meee stabliceh the initial time between reviews: 


2) and calculate the probability of a stock-out: 


RISK 


where: 


HEC 


n(r) 


meme Caleulate order cost = 


EEC 


ie 


stock-out cost per unit 

Gpder costyper order 

maximum inventory level 

holding cost per unit 

mean demand in an order lead time 
expected number of stock-outs per cycle 
price per unit 

annual demand 

review cost per review 

review period in years 


mean demand in an order cycle 


HD 
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4) Calculate the maximum inventory level using the 
initially determined probability ied Seeeck-our H(r). If 
this is a new iteration, compute a new H(r) based on the 
changed T value from step six. 

Consulting a Poisson CCDF for a mean demand in a 
lead time and an order cycle (M + D), find the E value 
corresponding to the point where the complementary cumulative 
Sieeene distribution (CCDF) is just greater than the RISK. 
This E value will be used to compute the rest of the cost 
ome tlOns. 


16 F 
Hor) = ¢f PCs) an lok 
Beate 


5) Find the expected number of stockouts per cycle from: 
mes) = ((M + D) * CCDF(E-1)) = (E * CCDF(E)) 


6) Compute the expected annual holding cost using the 


generated E from step two: 
molding cost = IC LE - > + )\M] 


7) Determine the expected annual stock-out costs: 


BE@ek-OutL Cost = kes 
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Byeeedd up the order, holding, and stock-out costs from 


step 3, 6, and 7 to get the total variable costs. 


AT 


mc = = +o0C a ele ee) 


Ki > 


Returning to step 3, a different T is chosen in order to 
calculate a new probability of stock-out. The T to be 
chosen will vary with each iteration as changes in total 
variable costs are observed. Finally, the correct maximum 
inventory level E and time between reviews T will be 
obtained when the iterations no longer result in a lower 
total variable cost. 

Since the purpose of this thesis is to compare the 
P-model with the Q-model (constrained to operate as the 


P-model) the costs of doing the review must be considered 


eventually. For the P-model the annual review costs are: 
Review cost = RE 


The total optimal costs calculated provide a baseline 
by which these costs can be compared with similar costs 
generated in the Q-model when that model is constrained to 
operate as the P-model. 

The major difficulty with the fixed order period model 
1s the system's limited responsiveness to a rapidly changing 
demand. Management has only limited knowledge, if any, of 


demand changes until the review period is reached and 


20 





shortages are more likely to result. Review costs are lower 
by having inventory examined at periodic intervals, but there 


is less information flow and control over inventory with this 


model than with the Q-model. 
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ee Leo UAT TONS 


A. BACKGROUND 

The simulations for the continuous review and fixed 
period models were developed to accurately represent the 
operation of an actual inventory system. Even though there 
are significant differences between the Q and P models, 
the major part of the simulations are common to both models. 
This chapter will describe the operation and assumptions of 
the simulations. The actual programs for the simulations 
are included as Appendix A and B. 

As mentioned in Chapter II, the Poisson distribution is 
used to generate demand. It is assumed that requisitions 
come in one at a time for a total quantity of one. A 
subroutine randomly generates the elapsed time between 
requisition arrivals (in days) based upon a mean expected 
daily demand. The times between arrivals are accumulated to 
provide a record of the total time of the simulation and how 
much time has passed since the last review. 

The time between reviews 1s used similarly in both models, 
but the time is input in different ways. In the continuous 
review model the time between reviews is read in from a data 
File so that this time can be altered to represent a fixed 
period constraint. The initial review time is always set at 


zero so that the system will be reviewed upon the receipt of 


Ze 








every requisition. This allows a comparison to be made 
between the theoretical optimal total cost--and the total 
cost of the simulation when the continuous review model is 
being operated in an optimal way. This check insures that 
the simulation is accurately representing the system. 
Succeeding times between reviews are then input based on 
how many days away from the optimal that the system is 


desired to be deviated. 


Eee COMPUTATIONS 

The fixed period model's initial time between reviews is 
calculated based on the initializing formulas outlined in 
Chapter II. This first review time will not necessarily 
represent the optimal time between reviews for the system. 
The optimal review time 1s computed using the iterative 
procedure described in Chapter II. This optimal review time 
1S input by the use of a predetermined percentage factor 
that sets the optimal time between reviews into the program. 
This percentage factor also allows the simulation of the 
P-model to be deviated away from the optimal for the purpose 
of examining alternative methods of operating the fixed 
period review system. 

The ability to vary the time between reviews is important 
because it is the basic way by which the two models can be 


compared. For example, the simulation can find the optimal 
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time between reviews in the P-model, and then fix that time 
into the Q-model. The resulting costs of both models can 
then be compared and analyzed. 

The key parameters in both models can also be altered in 
order to analyze different types of inventory items. A most 
important factor, mean annual demand, can be changed by using 
different values for the mean time between arrivals in the 
Poisson process (RLAM). For example, a time between arrivals 
of 0.5 (in days) would generate a mean number of 2 requisitions 
a day or 740 a year. A time between arrivals of 8.0 would 
cause a mean demand volume of about 45 per year. As the time 
between arrivals increases, and hence decreases the yearly 
mean demand, the variability of the Poisson process will also 
become more apparent. Other parameters such as holding, 
order, and shortage cost can also be changed. 

Reorders are generated differently by the models, but 
order quantities for the two systems are both based on 
inventory position. Inventory position is defined as a 
combination of the on-hand inventory level plus currently 
on-order minus any backorders. The continuous review system 
Will create a reorder when both the time between reviews 
has elapsed and the inventory position has fallen below the 
reorder point. The Q-model, when operated in a periodic 
review mode, has as a reorder quantity the difference between 
the maximum inventory position (sum of the reorder point and 


economic order quantity) and the current inventory position. 
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This order policy corresponds to what actually happens in 

the Navy's wholesale inventory system if the breach of a 
reorder point is not detected immediately due to some review 
time delay. Operating the system this way prevents the 
problem of inventory levels becoming negative and never 
rising above the reorder point again. A reorder in the fixed 
period review model is generated when the time between 
reviews has elapsed and the inventory position is less than 
the maximum established inventory level. 

Two subroutines are used to generate orders, keep track 
of when the orders are due to come in, and add to the on-hand 
inventory the appropriate quantity when an order has arrived. 
A constant lead time, which can be varied on different runs, 
is input at the beginning of the program. Any number of 
orders can be outstanding at a time. The time ramaining 
until each order will arrive is decreased, as days pass, 
until there is a simulated receipt. 

Theoretical costs for each simulation are computed as 
discussed in Chapter II. These theoretical costs are compared 
with the costs actually resulting from running the simulation. 
Simulation costs result strictly from what actually happens 
when the program 1s run under the inputed constraints and 
variables. For example, if five orders per year are 
generated then the total order cost will be five times the 


appropriate cost to order. 
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i eeCoUbAn toon OF CURVES 


A. BACKGROUND 

As originally stated, the purpose of this thesis is to 
find the "crossover" point in an inventory system which it 
1s advantageous to shift from a continuous review model to 
a periodic review model. The simulations, as outlined in 
Chapter III, provide the mechanism for finding this point by 
being able to run both the Q and P models under similar 
Memeare1Ons so that a comparison eee ewes et ing seOmad 
variable cost curves can be made. This chapter describes 
the details of the methodology by which these "crossover" 
points are found. 

There are a infinite number of possible combinations of 
holding, shortage, and order cost that could be used in the 
Simulations. Similarly, it is possible to go through endless 
variations of item price, lead time, and demand volume. As 
Mme@emonstraticon of the type of information that is desired 
from the simulations, an example follows in this chapter 
which will limit itself to one particular inventory item. 
The variable factors and costs that comprise this one 


inventory item are as follows: 


Item price: SO O0 

Holding rate percentage: w2ceper dollar per year 
Order cost: $175.00 per order 
Shortage cost: $350.00 per unit/year 
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The above parameters will remain constant throughout all 
Simulation runs. The only parameter that will be varied is 
the yearly number of requisitions (annual demand volume). 

A second parameter, lead time (in days), will also be varied 
but in a limited way. Three lead times will provide a 
sufficient number of curves to demonstrate the effects of a 
varying lead time on the crossover points of the two models. 
At this point the cost to do the review is ignored. Review 
cost impact on these curves will be analyzed in the next 


Simapter. 


Eee CURVE CONSTRUCTION 

It was necessary to run each simulation a number of times 
to get an average of the total annual variable costs for both 
models. In order to provide this "averaging" a minimum of ten 
Simulation runs were made with all factors held constant. If 
these runs showed low variability in the resulting cost data 
no further simulations were conducted. If high variability 
was encountered more simulations were run until a good 
consistent average of costs could be obtained. Furthermore, 
to provide an opportunity for the programs to produce an 
average of costs over a long period of time, the total length 
of each simulation run was fifty years. 

After examining data for several Navy managed inventory 
items it was observed that a rough average lead time for 


most stock numbers was in the range of five to seven quarters. 
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For this reason the first set of simulation runs used a 
lead time of 550 days. With this "set" lead time demand 
was selectively varied from two requisitions per year to 
one-hundred requisitions per year. 

Figure 1 is an example of the crossover point. Figure l 
is for a demand rate of twenty-five requisitions per year and 
shows the Q-model starting from its optimal total variable 
cost--the point where the inventory system 1s reviewed 
every day. As the Q-model is constrained to only review the 
system on a gradually increasing time basis, the total 
variable cost begins to increase. This increase in costs 
is primarily due to the higher number of stock-outs that 
are incurred because they are not detected until a review 
is made. Initially, the rise in costs are gradual, but as 
more and more days elapse between observations of inventory 
levels the costs rise at an increasing rate. 

In the P-model the costs of having a daily review are 
very high. This is due to the assumption that an item is 
consistently ordered up to a maximum inventory level every 
time the inventory falls below that specified level. 

Slowly, as the time between reviews increases, annual order 
costs go down until some optimal point where the total 
annual variable costs of the P-model are at a minimum. Even 
though order costs continue to decrease after this minimun, 


the further delay in reviews causes an increase in stock-out 
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costs. This increase in stock-out costs eventually 
outweighs the savings in order costs. 

The place where the variable cost curves intersect is 
termed the crossover point. In Figure 1 this point is a 
time between review periods of 99 days. If a inventory 
manager, uSing a continuous review system, allows the time 
between reviews to extend beyond 99 days then he/she is 
better off going to a fixed period review model to obtain 
a lower total variable cost. 

It 1S important to note that any time one of tnese 
models is not operated at its minimal cost point a penalty 
is being incurred by the inventory system. Even though the 
periodic review model in Figure 1 does not become preferred 
over the continuous review model until after 99 days, the 
continuous review system still is being operated non-optimally 
unless the system 1S reviewed every day. Conversely, for 
the P-model, using a time between reviews that is different 
from the bottom of the total variable cost curve at 160 days 
will result in a higher cost for the system--even though it 
is cheaper than using the Q-model. 

Figure 2 presents a set of crossover points for a lead 
time of 550 days and six different demand rates including that 
of Figure 1. It is easy to see that the crossover point 


decreases with increasing yearly demand volume. 
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From data such as provided by Figure 2 it is possible 
MeomeQhistruct a plot of crossover point times as a function 
Seecne mean annual demand. Figure 3 is such a plot for 
the information from Figure 2. For example, if demand is 
fifty requisitions per year it is favorable to switch 
inventory models if the time between inventory reviews slips 
beyond 62 days. Similarly, a low mean demand of only five 
per year would have a crossover point of over 202 days. 

Finally, other simulations were run in the same manner 
as outlined in this chapter but with changing lead times. 
One set of simulations was run with a lead time of three 
quarters (275 days) and another set with a lead time of 
mine quarters (825 days). Figure 4 presents these two new 
curves as well as the original curve of Figure 3. 

Figure 4 shows that as lead times decrease the total 
number of days to the crossover increases. For example, at 
a demand rate of fifty per year a lead time of 825 days has 
a crossover of 85 days, a lead time of 550 days has one of 
62 days, and a lead time of 275 days has one of 24 days. 
These results are quite logical since missing a reorder 
point would be expected to be more devastating in a long 
fead time situation. 

This chapter has presented a methodology for establishing 
the exact point that it is preferable to change inventory 


systems from continuous review to periodic review. It has 
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highlighted only one particular set of item cost data. 
However, this method does provide a potential way for an 
inventroy manager to decide if a change in inventory models 
is necessary if he/she finds themselves unable to review 


inventory levels upon the receipt of every requisition. 


So 





V. COST ANALYSIS 


A. INTRODUCTION 

As stated in Chapter IV, the simulation runs were limited 
to one inventory item. As such, any generalizations about 
variable cost behavior would only apply to inventory items 
with relatively similar cost parameters. However, there were 
trends observed with possible applicability over an entire 


inventory range. 


BeeeoFliIMAL MODEL COSTS 

In all cases the optimal total variable cost of the 
continuous review model (Qopt) was lower than the optimal 
variable cost of the periodic review model (Popt). Table I 
(columns 2 and 3) list these optimal costs in lead time 
categories by varying demand rates. The differences between 
these optimal costs (Col. 4) varied from 2 to 41% with the 
largest differences being experienced in the lower demand 
rates. The reasons for this trend is that P-model order 
costs, as a percentage of optimal total variable costs, 
increases as demand decreases. The Q-model order cost 
percentage remains relatively stable as demand changes. This 
causes the percentage difference in total variable costs for 
the two models to increase as demand falls. In addition, as 


lead time decreased the average differences between Qopt and 
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LEAD TIME AND DEMAND VARIABLE COST COMPARISONS 
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Popt consistently increased. The data indicated a linear 
relationship showing a 1.1% decrease in optimal cost 


differences for every 100 day increase in lead time days. 


Gree OSSOVER POINT COSTS 

Given constnat lead times the crossover point of the two 
models always increased as demand rates decreased. These 
periods between reviews were quite large, sometimes extending 
out to a year, for items with low demand. Conversely, highly 
demanded items experienced crossover points as short as 13 
Peer Lhis information is listed in column 5 of Table I and 
was used to plot Figure 4. 

The fact that a point exists where switching from a 
Q-model to a P-model will always result in a cost improvement 
does not alleviate the problem that both models are being 
operated in a non-optimal fashion. Table II highlights this 
by showing the Q and P model variable costs as a function of 
the review interval for just one item with a demand rate of 
25 and a lead time of 550 days. The crossover for this item 
meme 100 days. At this point the ratio between TVC for the 
Q-model and the TVC for the P-model is 1.00 and the total 
variable costs are both $1160. Operating the Q-model 
optimally at 0 days of review would generate only $901 in 
costs so allowing the period of review to slip to 100 days 
Will represent a 29% increase in variable costs to the 


Q-model system. Similarly, the P-model finds its optimal 
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SINGLE ITEM OPTIMAL VARIABLE COST COMPARISON 
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cost point at $943 (a review period of 160 days) representing 
a 22% loss in cost optimality at the crossover point. 

Table I continues with the analysis done above by listing 
the total variable costs at the crossover point, in column 6, 
for all simulated demand and lead times. In general, the 
higher the demand rate the greater the difference between the 
costs at the crossover point and the optimal costs of either 
model. This is significant when considering that the items 
with the higher demand are also the items with the highest 
variable costs. Straying from model optimality, in high 
demand situations, causes a larger cost penalty to the system 
than would be experienced by a lesser demanded item in the 


Soamencircumstances. 


eee LEW COSTS 

As mentioned earlier in this chapter, the optimal total 
variable costs of the Q-model appears to be always lower than 
the same cost for the P-model. However, there is a point 
where the increasing variable costs of the continuous 
review model (due to longer review periods) will be higher 
than the optimal variable cost of the periodic review model. 
In all cases this review interval was found to be much less 
than the crossover point. Table II demonstrates this is 
column 4 by comparing the ratio of changing Q-model variable 
costs to the constant P-model optimal cost. At a period of 


review of 25 days the ratio goes to 1.00 and it is there 
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that it 1s optimal from a cost standpoint to switch from the 
Q to the P-model. However, if this is done, the period of 
review will have to be 160 days, not 25 days, in order to 
imei benefit from the switch. 

All the variable costs in Chapter IV and Tables I and II 
do not consider the impact of the actual costs to do the 
review on the two models. For the data in Table II, the 
P-model will only by reviewed every 160 days. The Q-model, 
at its optimal point, is reviewed at a minimum upon the 
receipt of every requisition--or an average of 25 times a 
year. This imposes a greater cost to the continuous review 
system than is being experienced by the periodic review 
system. If it is true that, at a certain period of review, 
the P-model becomes cheaper to operate than the Q-model then 
there is an review cost that will made the optimal costs of 
Pemmemodels equal. 

For the data in Table II the review cost that will make 
both models equally optimal is $1.36. This $1.36 represents 
Meemort the total cost of the example inventory item ($300). 
If every time a manager examines inventory levels it costs 
the system greater than $1.36, then it is optimal from a 
cost standpoint to switch to the P-model and review every 
160 days. Conversely, if the review cost is less than $1.36, 
the Q-model is optimal with reviews done upon the receipt 


of every requisition. 
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Table I lists in column 7 the relevant review costs that 
will cause equal model optimality for varying demand and 
lead times. In general, items with demand greater than 25 
had more consistent, and lower, equality review costs than 
very low demanded items. 

The costs to do a review had no impact on the review 
interval at the crossover point. This is logical since, at 
the crossover point, both models are reviewed the equivilant 
number of times a year. The only change 1S an equal increase 


mmmGetal variable cost. 
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Veo UMMAR Vee) CONGEUS TONS 


A. SUMMARY 

The choice of an appropriate model to use in an inventory 
system will depend on a number of different factors. Some of 
these factors are related to intrinsic costs within the model 
itself. An example of these factors would be the order, 
Mewaame, and shortage costs that a particular activity might 
find applicable for the business conditions in which they 
operate. Other factors will also affect the choice of 
an inventory model. These factors include the actual size 
of the inventory system, the availability of personnel and 
computer resources, or the determination of what level of 
service was to be maintained. Regardless of which model is 
chosen there is an optimal way to operate that one inventory 
system. 

Sometimes it becomes necessary to deviate from the 
optimal operation of a particular inventory model. Ifa 
computer goes down it might not be possible to check stock 
levels upon the receipt of every requisition as is required 
by the "optimal" conditions of the continuous review model. 
If the mean demand increases dramatically the optimal set 
time between reviews in the periodic review model might be 


too long to adequately provide customer support. As these 
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various Omehtes ons change, a point can be reached where it 
is better to completely switch inventory systems rather than 
to continue distorting the present system. 

Two major models wien ee the basis by which most inventory 
systems operate. These models, the continuous review model 
and periodic review model, were analyzed in Chapter II. The 
purpose of this thesis was to find the point where the vari- 
able costs of these two models, under a given set of conditions, 
were equal. This area of equality was termined the "crossover" 
and it provided a point, for which, a switch of inventory 
models would be applicable from an annual variable cost aspect 
ignoring review costs. 

Simulations of each model were developed to find this 
crossover point. Chapter III describes the simulations and 
the underlying assumptions that make them work. Both simu- 
lations were built as elementary representations the United 
States Navy's wholesale inventory system under the conditions 
of the periodic and continuous review models. The programs 
are included as Appendices A and B. 

Using a given set of cost parameters a large number of 
simulations were run for both models at varying mean annual 
demand levels. A period of review constraint was imposed 
within the simulations to provide the change in variable 
costs that results from operating the models non-optimally. 


Once a steady state set of costs was obtained, the total 
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variable costs of the continuous and periodic review models 
were piloted as a function of the period of review for both 
models. The point where these cost curves intersect was 
the Mesesover point for that selected demand level. 

Using varying mean annual demand values, a series of 
crossover points was determined for a given set of other 
model parameters. The crossover point values were then 
plotted as a function of mean annual demand. 

The affect of lead time on the crossover point was also 
examined. Additional plots were developed from the simula- 
tion runs by changing lead time values. Chapter IV details 
the process by which the various curves were created. 
Chapter V then analyzed cost relationships between the two 
models for a given set of parameters. At the optimal point 
the Q-model total variable costs were always lower than 
optimal P-model variable costs. Operating either model at 
the crossover point was found to be inefficient and costly-- 
particularily for items with high demand. If review costs 
are imposed on each model, a dollar figure per review could 
be found that would equate the minimum variable costs of the 
continuous and periodic review models. Review costs had no 


impact on the review interval at the crossover point. 
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Bee CONCLUSIONS 

The decision of which inventory model is best for an 
activity must consider a number of variables in addition to 
those analyzed in this thesis. However, the procedures 
outlined above provide a mechanism for getting such informa- 
tion to the decision maker. The value of this information 
increases if the inventory system has to be operated 
non-optimally due to factors beyond the control of the 
inventory manager. Furthermore, if deviations from 
optimality do occur they must be addressed since these 
deviations both impose a penalty on the current system in 
higher variable costs and, at the same time, potentially 
offer opportunities for lowering costs by changing inventory 


Systems. 


Cees COMMENDATIONS 

All efforts should be made to operate inventory systems 
consistent with the optimality assumptions that created the 
inventory models. Any system constrained not to do so must 
be constantly analyzed to ensure that the correct inventory 
model is being used. A methodology, such as is presented in 
this thesis, should be available to the manager to assist 


in that decision process. 
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APPENDIX A 


CONTINUOUS REVIEW MODEL SIMULATION PROGRAM 


> 
ey 
= 
to 


CO.) CQ Oe GG O11 9 1 1 0) OO) OO) OOOO Oa as 


THES PROGRAM IS A SIMULATION OF A THEORETICAL INVENTORY 
SYSTEM WHERE A CONTINUOUS REVIEW MODEL IS CONSTRAINED TO 

Pee kENT TIMES OF FIXED PERIOD REVIEW. THE PROGRAM IS 
WRITTEN IN THE FORTRAN PROGRAMMING LANGUAGE USING THE WATFIV 
CeurPiLER 


Tit FOLLOWING IS A DESCRIPTION OF THE VARIABLES WITHIN THE 
PROGRAM. 


fee DEL F.G,H,1,0,K,1,0,R,S,1,U,V,W,X,Y,Z,AA,BB,VV,CC,GG, 
HH,SS, = ALL INTERNAL WORKING VARIABLES. 


AVGSOA = AVERAGE STOCK ON HAND AT REORDER 

BAKORD = TOTAL NUMBER OF BACKORDERS 

Seeex = DIGIT TO CAUSE NEW SETTIM VALUE TO BE READ 
DSEED = SEED FOR RANDOM NUMBER GrNERATOR 


EQOA = COMPUTED EOQ AFTER SIMULATION RUN 
Beeoly = ECONOMIC ORDER QUANTITY COMPUTED BEFORE SIMULATION 
RUN 


Pee bo) = EXPECTED NUMBER OF BACKORDERS THEORETICAL 
eieeol = HOLDING COST PERK ITEM PER YEAR 

B@ieek = HOLDING RATE PERCENTAGE (0<X<1) 

Peeeeie = TIME BETWEEN ARRIVALS TAKEN FROM ARRAY IR 


fee UMBER OF REQUISITIONS 

INVLEV = INVENTORY LEVEL ACTUALLY ON-HAND 

INVPOS = INVENTORY POSITION (ON-HAND + ON-ORDER) 

IR = ARRAY OF POISSON GENERATED TIME BETWEEN ARRIVALS 
LEADTM = MEAN ORDER LEADTIME IN DAYS 

MAXA = COMPUTED MAXIMUM INVENTORY AFTER SIMULATION RUN 
MAXINV = MAXIMUM INVENTORY LEVEL 


Niee= NUMBER OF ITEMS IN ARRAY IR 
Bee akRAY THAT KEEPS TRACK OF THE NUMBER OF REQUISITIONS 
Memes okRAY THAT KEEPS TRACK OF LEADTIMES FOR ORDERS 


ONORD = TOTAL AMOUNT ON ORDER AT ANY ONE TIME 
MemerakkAY THAT KEEPS TRACK OF ON ORDER QUANTITIES 
ORDCNT = TOTAL NUMBER OF ORDERS DURING THE SIMULATION 
Oleeol = ORDERING COST PER ORDER 

ORDYR = AVERAGE ORDERS PER YEAR 

RLAM = MZ&AN ARRIVAL RATE IN DAYS 

PeolykR = AVERAGE REQUISITIONS PER YEAR 

REORD = REORDER POINT 

REVCNT = COUNTS TOTAL NUMBER OF REVIEWS 

REVCST = COST PER REVIEW 

REVYR = AVERAGE REVIEWS PER YEAR 

PeviiM = CUMULATES TIME BETWEEN REVIEWS IN DAYS 
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Pepolk = SAFETY STOCK THEORETICAL 

SETTIM = THIS 15 TOTAL TIME ALLOWED BETWEEN REVIEWS 
Saoeol = SHORTAGE COST PER SHORTAGE 

SIMRUN = TOTAL LENGTH OF SIMULATION RUN IN DAYS 

SMA = SUPPLY MANAGEMENT AVAILABILITY 

STKOA = TOTAL STOCK ON HAND AT ORDER ARRIVAL 

pexwk = AVERAGE STOCK-OUTS PER YEAR 

eee ey = TOTAL LEADTIMES SUMMED FOR ALL ORDERS 

SUMORD = TOTAL NUMBER OF ORDERS 

TOTA = TOTAL YEARLY COST OF RUN WITH REVIEW CONSTRAINT 
meaeo = TOTAL OPTIMAL COST BASED ON INITIAL INPUT VALUES 
meme = KEEPS TRACK OF REAL TIME OF SIMULATION 

ieeoo = TOTAL NUMBER OF STOCK-OUTS ON RUN 

TVCO = TOTAL VARIABLE COST OPTIMAL 

TVCR = TOTAL VARIABLE COST SIMULATION RUN 

maps = FACTOR TO COMPUTE YEARLY COSTS 

DECLARE VARIABLES AND ARRAYS 


REAL TVCR,RVCO,REQNYR,REVYR,ORDYR,SOYR, YEAR, SAFSTK , EXPBO 
REAL HOLCST ,PRICE,HOLPER,A,FLAM,ORDCST , SHOCST,B,EOQQTY 
REAL MAXINV,C,BAKORD,REORD,EOQA,MAXA,D,E,SMA,REVCST 

REAL TOTCO,Z,Y,X,W,V,U,TOTA,AA,BB,H,G,R,S,T,CC,GG,HH,SS 
REAL STKOA,AVGSOA 


miele ONORD,J ,K,£,SUMORD,SUMLD,NR,IER,SIMRUN ,TOTIME 
mieonk HITIME ,SETTIM,REVIIM, INVLEV ,TOTSO,INVPOS,ONORD 
ieeonk REVCNT,ORDCNT ,LEADTIM,CHECK,N 

EGER IR€2500) ,OC(100) ,0Q(100) ,OL(100) 


POUBLE PRECISION DSEED 


INITIALIZE VARIABLES AND SET COST AND RUNTIME FIGURES. 
READ THE PERIOD OF FIXED REVIEW AND MAKE A CHECK FOR 
END OF DATA RUN 

READ(5 ,810)SETTIM, CHECK 

DSEED = 35110.0D0 

HOLCST = 0.0 

EOQQTY = 0.0 

MAXINV = 0.0 

oe) 0.0 

Ee 0.0 

om-) 0.0 

REORD = 77.0 

EXPBO = .034 

BAKORD = 0.0 

NR = 2500 


BeecCE = 100.0 
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CDEC O17. 


C1 7G): 


HOLPER 
ORDCST 
SmoCST 
RLAM = 
SIMRUN 
LEADTM 
eee IME 
BETiME 
REVTIM 
SIMCNT 
REVCNT 
ORDCNT 
SUMO RD 
SUMLD 

ONORD 

mo S0 

STKOA 


29 
iS. 0 
aU). U 


14600 


Hott th tt th ht ae 
SO © Ok & 


SS Ge 


0 
Meno ORDER ARRAYS 


DO 5 J=1,100 
Gear) =0 
0Q(J)=0 
OL(J)=0 

CONTINUE 


BeGlN PREPROGRAM CALCULATIONS WITH THE CALCULATION 
Meee ECONOMIC ORDER QUANTITY (E0Q). DEMAND IN THIS 
Moe ts FEGURED BY DIVIDING DAYS IN THE YEAR BY THE 
MEAN ARRIVAL RATE (365.0/RLAM). 


HOLCST = (PRICE * HOLPER) 

A = 2.0 * (365.0/RLAM) 

B = (SHOCST * EXPBO) + ORDCST 
@em- <A * B)/HOLCST 

EOQQTY = SQRT(CC) 


eee AX EMUM INVENTORY LEVEL IS COMPUTED BY ADDING THE EOQ 
QUANTITY TO THE REORDER POINT. THE INITIAL INVENTORY LEVEL 
Pee VENTORY POSITION IS SET AT HALF THE EOQ QUANTITY PLUS 
elmewromDER POINT. SAFETY STOCK IS COMPUTED BY 

PPE RACTING DEMAND IN A LEADTIME FROM THE REORDER POINT. 


meee oORTCCSHOCST) / CHOLCST + SHOCST)) 


MAXINV = EOQQTY + REORD 

MiGeavy = LEFIXCREORD) + IFIXCEOQQTY)/2 
INVPOS = INVLEV 

SAFSTK = REORD - (FLOAT(CLEADTM) /RLAM) 
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CO) GO 


350 INVLEV 


QO Bo 


(Go <)) C2 


Ore @ @ (oGa7eo 


mem ooON ARRIVALS ARE GENERATED USING THE NAVAL POSTGRADUATE 
PenoOl Ss IMSL SUBROUTPNE GGPOS. THIS SUBROUTINE MUST BE 
imyOKED IN IBM 3033 CMS BY USING THE STATEMENT 

PebesAL MACLIB IMSLDP™ PRIOR TO THE RUNNING OF THE PROGRAM. 


CALL GGPOS(RLAM,DSEED ,NR,IR,IER) 
BEGIN SIMULATION 
DO 50 I=1,2500 


CHECK TO SEE IF SIMULATION RUN TIME Pome vrn CUNEATE =f IME 
BETWEEN REVIEWS 


PE ChLORME. GE. olMRUN) GO T0500 
Feels blll IR(T) 

ASO nMeie iO ee eet iE 

Ry eh Rev Met he EME 


Mee, TO SEE IF SYSTEM £S OUT OF STOCK AND RECALCULATE 
Mey ENTORY LEVELS BASED ON NEW REQUISITION. 


EPPCENVEEY -GT.0) GOMmO 300 

BOM SOmae LOLS O +1. 
INVLEV - l 
INVLEV + ONORD 


INVPOS 


Seren FO SEE IF TIME BETWEEN KEVIEWS IS FINISHED 
Pee Reobl THE TIME BETWEEN REVIEWS. 


PE CSELLOIM. GE Rev Er) GCOemer 400 
REVTIM = 0 
REVCNT = REVCNT + 1 


@meek TO SEE IF INVENTORY POSITION IS BELOW THE 
Rete ORDER POINT. IF 50, CALL UP REORDER SUBROUTINE. 


PECtivE@s.eh. LEexCREORD)) GO TO 400 
CRDeCNE = ORDCNE + 


GENERATE A REORDER 
CALL ORDER(MAXINV,ONNORD,LEADTM, INVLEV,0C,0L,0Q,INVPOS) 


PemORE GENERATING THE NEXT REQUISITION CALL UPDATE SUBROUTINE 
Ponce WiLL UPDATE INVENTORY LEVELS AND ONORDER LEADTIMES 


90 








400 


on 


SOV OG Seo a oe 


on 
(om 
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QO ee 


CALL UPDATECONORD,HITIME,SUMORD,INVLEV,SUMLD,0C,0L,0Q, 
INVPOS ,ORDCNT , STKOA ,AVGSOA) 


SLART A NEW REQUISITION INTO THE SYSTEM 
CONTINUE 
Seer eh FINAL STATISTICS AS YOU EXLT PROGRAM 


meee ec. TON COMPUTES THE ACTUAL EOQ QUANTITY, MAXIMUM 
ENVENTORY LEVEL, AND SUPPLY MANAGEMENT AVAILABILITY THAT 
MoUGED HAVE RESULTED FROM THIS SIMULATION 


YEAR = 365.0/FLOAT(SIMRUN ) 
ORDYR = FLOATCORDCNT) * YEAR 
D = YEAR * FLOAT(T) 

fe FLOATCTOTSO) * YEAR 

eee 0 / ORDYR ) 

fie = (SHOCST * GG) + ORDCST 
eree- (2.0 * DPD * HH) /HOLCST 


EOQA = SQRT(SS) 
MAXA = EOQA + REORD 


SMA = (FLOAT(I) - FLOATCTOTSO) )/FLOAT(I) 


MieosoeCliON COMPUTES THE TOTAL THEORETICAL COSTS OF THIS 
Semel BASED SOLELY ON THE INITIAL INPUT VALUES. NOTE THAT 
Pere COSTS WILL BE INCURRED ON EVERY REQUISITION 


(RLAM * 365.0) * PRICE 

(ORDCST * (365.0/RLAM)) / EOQOQOTY 

HOLCST * ((EOQQTY/2.0) + SAFSTK + EXPBO) 
(365.0/RLAM) / EOQQTY 

SHOCST * EXPBO 

REVCST * (365.0/RLAM) 

W=WwWe#y 

TOTCO = Z+Y+X+U+VV 

TVCO = Y+X+VV 


THIS COMPUTES THE TOTAL YEARLY COSTS OF OPERATING THIS 
feel baoky ON THE ACTUAL RESULTS OF THE SIMULATION. 


PRICE * (365.0/FLOAT(SIMRUM)) * FLOAT(TI) 
(ORDCST*FLOAT (ORDCNT) ) *(365.0/FLOAT(SIMRUN) ) 
REORD - (FLOAT(LEADTM)/365.0) * (FLOAT(I) * YEAR) 
(EOQA/2.0) + R+E 

HOLCST * Q 
= FLOAT (TOTSO) *SHOCST*(365.0/FLOAT(SIMRUN) 
= (FLOAT(REVCNT)*REVCST)*(365.0/ FLOAT (SIMRUM) ) 
TOTA = T+S+G+AA+BB 
S+G+AA 
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meboeoechlON COMPUTES SUMMARY STATISTICS FOR THE 


SIMULATION RUN 


REQNYR = FLOAT( 


ORDYR = FLOATCORDCNT) * 


I) * YEAR 


pein 


SOYR = FLOAT(TOTSO) * YEAR 
REVYR = FLOAT(CREVCNT) * 


YEAR 


eon olALTEMENTS WRITE THE SUMMARY STATISTICS 


WRITE(6 ,600) 


WRITE(6, 700 )EQQQTY ,MAXINV,Z,Y,X,VV,U,TOTCO,TVCO,AVGSOAV 


WRITE(6 ,800) 


WRITE(6,380)SETTIM 


Pele (6, 900)EQOQA,MAXA,T,S,G,AA,BB,TIOTA,TVCR,SMA 


WRITE(6,300) 


mie 6,200) SAITIME ,TOTIME,REVTIM, REVCNT, INVLEV, 
INVPOS, ORDCNT, SUMORD, ONORD, SUMLD, TOTSO 


Pte (6,305) 


eee (653/75) REQNYR,REVYR,ORDYR,SOYR, SETTIM, eee, 
eC Eye hin Co Bee ORDCST, CHOCST, RLAM, SIMRUN 


moms oECIION FORMATS THE READ AND WRITE STATEMENTS 


BORMATCIX,"RUN WITH TIME BETWEEN REVIEWS SET AT: 


BORMATCL217,/) 


FORMAT (2X,'REQN! 


i Rey Ci six. 


1X,'SUMORD'1X, 
FORMAT(1X, "THEORETICAL"! 


iPkeOsT ' 
poled @S key 


ylx,'HITIME',1X,'TOTIME! 
INVLEV'! ,1X,'INVPO 
'ONORD' ,1X,'SUMLD! 
ie nOCOT wT. 5x. 
5x, 1ORDGST UM eu HOmesT 
(So Ores men ere 


FORMAT(12X,10F11.2,/) 
FORMAT (1X, 'SIMULATION 


'ITCOST', 
5X,'REVCST! 


' 2X, 'EOQOTY' 55x, 
'HOLCST! 
So ee nC 


5X,'ORDCST' 55x, 
5 erect 


memiverCi2xX,9F11.2,Fil.4,/) 


mee MAT (215) 


FORMAT (5X,'REQNYR!' ,5X,' 


6X,'SET! 
'ORDCST! 


FORMAT(1X,'MAC INVENTORY: 
Bee. 3X. 
F8.4,3X, 


Fe ei 


, 7X, 'LEADTM' ,5X, 
,5X,'SHOCST' 55x," 
FORMAT (4F11.2,2111,5F11.2,I111,//) 


eyavole hyd 


SLOCK: 


"GENERATED B/O: 


fern Ce 


,F8.2,3X, 


f 
3 


RLAM ! 


He!) 


,1X,'REVTIM', 
S',1X,'ORDCNT 


I aeeorso)) 


"MAXINV! ,5X, 


'MAXINV! ,5X, 


,5X,'SHOCST', 


'5X,' SMA 


,5X,' ORDYR',5X,' SOYR 
,5X,'REVCST' 5X, 
,5X,'SIMRUN') 


"REORDER POINT” ' 
SHORI26 3 


t 
3 


,5X,'SHOCST' 55x, 
,5X, 'AVGSOA') 


oe 


t 


3 


eee ele) BAO: 


READ IN A NEW TIME BETWEEN REVIEWS AND CHECK DIGIT 


ees, oLO)SETTIM CHECK 


mecewpCK.Gl.0) 


COM LOM ns 


a2 


3 


PNOrs/57 2 


, 


3 
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60 


65 


68 


i Gana 6 2 CD 0D 


80 


70 


BoeeD = DSSED + 112.0D0 


2 LOP 

END 

END OF PROGRAM 

BEGINNING OF SUBROUTINE SECTION 

ees SUBROUTINE COMPUTES A NEW ORDER. THE QUANTITY ORDER- 


Pemeiol; BE THE DIFFERENCE BETWEEN THE INVENTORY POSITION 
AND THE MAXIMUM INVENTORY LEVEL 


SUBROUTINE ORDER(MAXINV,ONORD,LEADTM, INVLEV,0C,0L,0Q,INVPOS) 
INTEGER 0C(100) ,OL(100) ,0Q(100) 
INTEGER INVLEV,ONORD,LEADTM,INVPOS 
REAL MAXINV 
DO 60 J=1,50 
IF(OC(J).EQ.0) GO TO 65 


CONTINUE 
GO TO 68 
a 
OC(K)=1 
O0Q(K)=(IFIX)MAXINV-INVLEV-ONORD 
. OL(K) =LEADTM 
RETURN 
END 
THIS ROUTINE UPDATES LEADTIMES AND ADDS TO INVENTORY LEVEL 
AN ONORDER QUANTITY IF IT IS TIME FOR THE ORDER TO COME IN. 
IT ALSO SUMS THE REQUISITIONS OUTSTANDING, THEIR QUANTITIES 
AND LEADTIMES 
SUBROUTINE UPDT(ONORD,HITIME,SUMORD,INVLEV,SUMLD,0C,OL, 
0Q,INVPOS ,ORDCNT,STKOA,AVGSOA) 


INTEGER ONORD,HITIME,SUMORD,INVLEV,SUMLD 
INTEGER INVPOS,N,ORDCNT 
INTEGER 0C(100) ,OL(100) ,0Q(100) 
REAL STKOA,AVGSOA 
DO 80 J=50 
TE OCG EO mOn@p Ol) )-HITIME 
CONTINUE 
N=51 


eou70) J=1. 50 


TE COCQW)  Giier@mebcd ),LE.HITIME) N=J 
CONTINUE 
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RETURN 
END 


IF(OQ(N).LE.0) GO TO 73 


ONORD=0 
SUMQTY=0 


SUMLDqQO 


STKOA = STKOA + FLOATCINVLEV) 
AVGSOA = STKOA/FLOAT(CORDCNT ) 
INVLEV=INVLEV+0Q(N) 

OC(N)=0 

OQ(N)=0 

OL(N)=0 


DO 90 J=1,50 


CONTINUE 


K=J 

ONO RD=ONORD+0Q(K) 
SUMO RD=SUMORD+0C(K) 
SUMLD=SUMLD+t+OL(K) 


INVPOS=INVLEV+ONORD 


ieee OLLOWING EXAMPLE DATA ENTRY IS THE TIME BETWEEN 
Pew AND A CHECK DIGIT FOR END OF RUN, TO STOP 
mEWLING VALUES ARE A CHECK DIGIT OF ZERO 


ENTRY 
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FIXED PERIOD REVIEW MODEL SIMULATION PROGRAM 


e 
S&S 
WO 


meso FROGRAM IS A SIMULATION OF A THEORETICAL INVENTORY 
SrolEM WHERE A FIXED PERIOD REVIEW MODEL IS USED BASED 
mien OPTIMAL TIME BETWEEN REVIEWS. THIS OPTIMAL TIME 
BETWEEN REVIEWS CAN BE OPTIONALLY CHANGED BY A PERCENTAGE 
Pienor. TO DETERMINE THE CORSTS OF OPERATING THE SYSTEM 

IN A NON-OPTIMAL MANNER. THE PROGRAM IS WRITTEN IN THE 
FORTRAN PROGRAMMING LANGUAGE USING THE WATFIV COMPILER 


THE FOLLOWING IS A DESCRIPTION OF THE VARIABLES WITHIN THE 
PROGRAM 


Peewee) EF Gst.1,0.K,1,0,R.$,1,U,V¥.W,X,Y.2,4A,BB,VV,CC,GG, 
HH,SS, = ALL INTERNAL WORKING VARIABLES. 


AVGREQ = AVERAGE REORDER POINT DURING SIMULATION RUN 
AVGSOA = AVERAGE STOCK ON HAND AT ORDER ARRIVAL 
hewORD = TOTAL NUMBER OF BACKORDERS. 

CHANGE = % CHANGE OF REVIEW PERIOD FROM TOPT 

mmeew = DIGIT TO CAUSE NEW SETTIM VALUE TO BE READ 
freee, = UNITS DEMANDED IN A REVIEW CYCLE 


Mee = YEARLY THEORETICAL DEMAND 

DSEED SEED FOR RANDOM NUMBER GENERATOR 

meee = COMPUTED EOQ AFTER SIMULATION RUN 

meeely = ECONOMIC ORDER QUANTITY COMPUTER BEFORE SIMULATION 
RUN 


ieeebo = EXPECTED NUMBER OF BACKORDERS THOERETICAL 
meeeol = HOLDING COST PER ITEM PER YEAR 
HOLPER = HOLDING RATE PERCENTAGE (0<X<1) 
PME = TIME BETWEEN ARRIVALS TAKEN FROM ARRAY IR 


IT = NUMBER OF REQUISITIONS 

INVLEV = INVENTORY LEVEL ACTUALLY ON-HAND 

INVPOS = INVENTORY POSITION (ON-HAND + ON-ORDER) 

IR = ARRAY OF POISSON GENERATED TIME BETWEEN ARRIVALS 


mewbey = LEADTIME OF REVIEW CYCLE IN DAYS 
meee = UNITS DEMANDED IN AN ORDER LEADTIME 
EEADTM = MEAN ORDER LEADTIME IN DAYS 


ieee = COMPUTED MAXIMUM INVENTORY AFTER SIMULATION RUN 
MAXINV = MAXIMUM INVENTORY LEVEL 


few= NUMBER OF ITEMS IN ARRAY IK 
Wiee= ARRAY THAT KEEPS TRACK OF THE NUMBER OF REQUISITIONS 
imeem ak RAY THAT KEEPS TRACK OF LEADTIMES FOR ORDERS 


eNekRD = TOTAL AMOUNT ON ORDER AT ANY ONE TIME 

woe —- ARRAY THAT KEEPS TRACK OF ON ORDER QUANTITIES 
ORDCNT = TOTAL NUMBER OF ORDERS DURING THE SIMULATION 
@ieeolr = OKDERING COST PER ORDER 


So 
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ORDYR = AVERAGE ORDERS PER YEAR 

PRICE = ITEM PRICE PER UNIT 

RLAM = MEAN ARRIVAL RATE IN DAYS 
REQNYR = AVERAGE REQUISITIONS PER YEAR 
REORD = REORDER POINT 


REVCNT = COUNTS TOTAL NUMBER OF REVIEWS 

Peer eol = COST PER REVIEW 

heey rR = AVERAGE REVIEWS PER YEAR 

fmevt iM = CUMULATES TIME BETWEEN REVIEWS IN DAYS 
BereolK = SAFETY STOCK THEORETICAL 

Peer = THIS IS TOTAL TIME ALLOWED BETWEEN REVIEWS 
SHOCST = SHORTAGE COST PER SHORTAGE 

SIMRUN = TOTAL LENGTH OF SIMULATION RUN IN DAYS 


Sie = SUPPLY MANAGEMENT AVAILABILITY 
seek = AVERAGE STOCK-OUTS PER YEAR 
Seite = TOTAL LEADTIMES SUMMED FOR ALL ORDERS 


eelionD = TOTAL NUMBER OF ORDERS 
TOLEAD = LEADTIME DAYS - REVIEW AND ORDER CYCLE 
MeevyyR = THEORETICAL REVIEWS PER YEAR 


mer = OPTIMAL TIME BETWEEN REVIEWS IN A YEAR 
imeETe = OPTIMAL TIME BETWEEN REVIEWS CHANGED BY % FACTOR 
mete = TOTAL YEARLY COST OF RUN WITH REVIEW CONSTRAINT 
moECO = TOTAL OPTIMAL COST BASED ON INITIAL INPUT VALUES 
Seeeeir, = KEEPS TRACK OF REAL TIME OF SIMULATION 
miso = TOTAL NUMBER OF STOCK-OUTS ON RUN 
memes = TOTAL UNITS DEMANDED IN A LEADTIME 

eee bon es Le Or tials 
TVCR = TOTAL VARIABLE COST SIMULATION RUN 

ee eon LO COMPUTE YEAKDY COstis 


DELCARE VARIABLES AND ARRAYS 


REAL TVCR,TVCO,REQNYR,REVYR,ORDYR,SOYR,YEAR,SAFSTK,EXPBO 
REAL HOLCST,PRICE,HOLPER,A,RLAM,ORDCST ,SHOCST,B,EOQQTY 
REAL MAXINV ,C,BAKORD,REORD,EOQA,MAXA,D,E,SMA,REVCST 

Pe TOTCO,Z,Y,X,W,V.,U,TOTA,AA,BB,H,G,R,$,1T,CC,GG,HH,SS 
REAL DEMAND ,TOPT ,TOPTC ,THRVYR,CYCDM,LEADDM,TOTUTS , AVGREQ 
REAL STKOA,AVGSOA 


INTEGER ONORD,J,K,1I,SUMORD,SUMLD,NR,IER,SIMRUN , TOTIME 
INTEGER HITIME,SETTIM,REVTIM,INVLEV ,TOTSO,INVPOS,CHECK 
INTEGER REVCNT,ORDCNT , LEADTM, LEADCY,TOLEAD,N 

INTEGER IR(2500) ,0C(100) ,0Q(100) ,OL(100) 


BOwBLE PRECISION DSEED 


PPI TALIZE VARIABLES AND SET COST AND RUNTIME FIGURES. 
mEeeeeCiED BACKORDERS AND MAXIMUM INVENTORY ARE DETERMINED 
PeesRATELY USING AN ITERATIVE PROCEDURE BASED ON THE 
FeeeeBiLitTyYy OF A STOCK-OUT, AND DEMAND IN A CYCLE AND 
Peet iME WETHIN A POISSON DISTRIBUTION 
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Cra) @ 


Om 2 


1 @) @  (ea 


DSEED 79984 .0D0 
CHANGE = 1.0 

MAXINV = 52.0 

EXPBO eZ 

DSEED 29146 .0 

NR = 2500 

meer = 100.0 

mor ekR = .23 

mem@sl = HOLPER * PRICE 
Oem est = 175.0 
eeewol = 390.0 
eevCol = 9.00 

RLAM = 1.0 

SIMRUN = 3650 
DEMAND = 365.0/RLAM 
LEADTM = 82 


THE FOLLOWING DETERMINES THE OPTIMAL TIME BETWEEN REVIEWS, 
eee eVIEWS PER YEAR, AND THE DAYS IN A CYCLE. 


MET = SORTC((2.0 * ORDCST)/( DEMAND * HOLCST) ) 
meene = TOPT * CHANGE 
MaayYR = 1.0/TOPTC 

MeveDCY = IFIX(TOPIC * 365.0) 
pope = LFIXCTOPTC * 365.0) 
REORD = 0.0 

ROLIME = 0 

REVTIM = 0 

REVCNT = 0 

ORDCNT = 0 

SUMORD = 0 

SUMLD = 0 

ONORD = QO 

OoTKOA = 0.0 

TOTSO = 0 


ZERO! ORDER ARRAYS 


DD 5 J=1,100 
OC(J)=0 
OQ(J)=0 
oni) =o 

CONTINUE 


BEGIN PREPROGRAM CALCULATIONS WITH THE CALCULATION 
meets, ECONOMIC ORDER QUANTITY (CEOQ). DEMAND IN THIS 
Woe Los PLGURED BY DIVIDING DAYS IN THE YEAR BY THE 
MEAN ARRIVAL RATE (365.0/RLAM). 


a) 





COO) Oe eo 


OCs" @ 


eer @) re) 


Bey G62 


SoU 


7: 


A = 2.0 * (365.0/RLAM) 

B = (SHOCST * EXPBO).+ ORDCST 
CC = (A * B)/HOLCST 

EOQQTY = SQRT(CCO 


ie NUMBER OF UNITS DEMANDED IN A CYCLE AND LEADTIME 
PoeecONPUIED HERE. INVENTORY LEVEL IS SET AT THE MAXIMUM 
MPENTORY LEVEL. SAFETY STOCK IS DETERMINED BY SUBTRACTING 
DEMAND IN A CYCLE AND A LEADTIME FROM THE MAXIMUM INVENTORY 
evel. 


CYCDM = FLOAT(CLEADCY)/RLAM 


LEADDM = FLOAT(LEADTM) /RLAM 
ments = CYCDM + LEADDM 
INVLEV = IFIX(MAXINV) 
INVPOS = INVLEV 

SAFSTK = MAXINV - TOTUTS 


POISSON ARRIVALS ARE GENERATED USING THE NAVAL POSTGRADUATE 
eeweol'S IMSL SUBROUTINE GGPOS. THIS SUBROUTINE MUST BE 
mrO.e> EN IBM 3033 CMS BY USING THE STATEMENT 

BemoBAl MACLIB IMSLDP" PRIOR TO THE RUNNING OF THE PROGRAM. 


CALL GGPOS(RLAM,DSEED,NR,IR,IER) 
BEGIN SIMULATION 
Comso 1=1,2500 
CHECK TO SEE IF SIMULATION RUN TIME IS OVER, CUMULATE TIME 


BETWEEN REVIEWS 
br CrORIMESGE oi RUNS GO TO 500 


EGE EME eee le) 
OT Eee kO leet ae EME; 
Fey eee le tol eee Lt ME: 


eee TO SEE IF SYSTEM ES OUT OF STOCK AND RECALCULATE 
BNVENTORY LEVELS BASED ON NEW REQUISITION . 


LPCINVEEV Gi Ops GOnLo ss 50 
TOTSO ="10850 * i 

INVLEV = INVLEV - 1 

INVPOS = INVLEV + ONORD 


Meeewlowoctk LF TIME BETWEEN REVIEWS IS FINISHED 
Pero RESEL THE TIME BETWEEN REVIEWS AND ORDER TO 
MAXIMUM INVENTORY LEVEL. DETERMINE THE AVERAGE 
Mev OF ENVENTORY AT THE ORDER. 
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= 
© 


on 
CO (2S GG) C.O Gree @ Oa) @ "(3 -@) 


on 
© 


OO: .@ 


Pe Combi ck. REVEEMOSGO TO 400 


REVTIM = 0 
REVCNT = REVCNT + 1 
ORDCNT = ORDCNT + 1 


REORD = REORD + FLOAT(INVLEV) 
AVGREQ = REORD/FLOAT(REVCNT) 


GENERATE A NEW ORDER 
ek, ORDERCMAXINV ,ONORD, LEADTM , INVLEV,0C,0L,0Q) 


Peeks GENERATING THE NEXT REQUISITION CALL UPDATE SUBROUTINE 
aot WELL UPDATE INVENTORY LEVELS AND ONORDER LEADTIMES 


eon UPDATE CONORD,HITIME,SUMORD,INVLEV ,SUMLD,0OC,O0OL,0Q, 
INVPOS ,;ORDCNT ,STKOA,AVGSOA) 


peek! A NEW REQUISITION INTO THE SYSTEM 
CONTINUE 
Seeeere FINAL STATISTICS AS YOU EXIT PROGRAM 


moreso otCliON COMPUTES THE ACTUAL EOQ QUANTITY, MAXIMUM 
MmewlORY LEVEL, AND SUPPLY MANAGEMENT AVAILABILITY THAT 
Pew HAVE RESULTED FROM THIS SIMULATION 


fee = 365.0/FLOAT(SIMRUN) 
SlpyR = FLOATCORDCNT) * YEAR 
De YEAR * FLOAT(TI) 

eee | LOATCTOTSO) * YEAR 

meee «6( 1 .0/ORDYR) 

fee (SHOCST * GG) + ORDCST 


SS = (2.0 * D * HH)/HOLCST 
EOQA = SQRT(SS) 
MAXA = MAXINV 


SMA = (FLOAT(I) - FLOAT (TOTSO)) /FLOAT(T) 


Bee seOr TON COMPUTES THE TOTAL THEORETICAL COSTS OF THIS 
eeaetbM BASED SOLELY ON THE INITIAL INPUT VALUES. NOTE THAT 
ioewW COSTS WILL BE INCURRED ON EVERY REQUISITION 


(© 55.0/RLAM) = PRICE 

end, LOrlC) OR DCsT 

Sever Alby sTOrle)/ 2.0) te 0DEMaNb. =. (FLOAT(LEADTM)/365.0)) 
HOLCST * (MAXINV - W) 

Gio, NOPTC) * SHOCST “ EXPEG 

ere POrPTC) “« REVCST 

meme = ZryYt+xtUt+Vv 

mye = YtxXtV 


OS 





Ca & 


© CGO) 


OOO 


Crea 


700 





toe osChlON COMPULES THE TOTAL YEARLY COSTS OF OPERATING 
ieomorolEM BASED ON THE ACTUAL RESULTS OF THE SIMULATION. 


T = PRICE * (365.0/FLOAT(SIMRUN) * FLOAT(TI) 

S = (ORDCST*FLOAT(ORDCNT) )*(365.0/FLOAT (SIMRUN) ) 
R = (YEAR * FLOAT(I) * TOPTC)/2 

Q = (YEAR * FLOAT(I)) * (FLOAT(LEADTM)/365.0) 

G = HOLCST * (MAXINV - (R + Q)) 

AA = FLOAT(TOTSO)*SHOCST*(365.0/FLOAT(SIMRUN) ) 
BB = (FLOAT(REVCNT)*REVCST)*(365.0/FLOAT(SIMRUN) ) 
TOTA = T+S+G+AA+BB 

TVCR = S+G+AA 

YEAR = 365.0/FLOAT(SIMRUN) 

THIS SECTION COMPUTES SUMMARY STATISTICS FOR THE 
SIMULATION RUN , 


REQNYR = FLOAT(I) * YEAR 
SOYR = FLOAT(TOTSO) * YEAR 
REVYR = FLOAT(REVCNT) * YEAR 


iteor SlATEMENTS WRITE THE SUMMARY STATISTICS 


WRITE(6,380)SETTIM 

WRITE(6,600) 

WRITE(6,700)EOQQTY ,MAXINV,Z,Y,X,V,Y,TOTCO,TVCO,AVGSOA 

WRITE(6, 800) 

WRITE(6,900)EQQA,MAXA,T,S,G,4A,BB,TOTA,TVCR, SMA 

WRITE(6,300) 

WRITE(6,250)1,HITIME ,TOTIME,REVTIM,REVCNT,INVLEV, 
INVPOS ,ORDCNT , SUMORD ,ONORD,SUMLD,TOTSO 

WRITE (6 ,255)MAXINV ,AVGREQ, SAFSTK , EXPBO ,GG 

WRITE(6,310)CYCDM,LEADDM,THRVYR 

WRITE(6,305) 

WRITE(6,375)REQNRY,REVYR,ORDYR,SOYR, SETTIM,LEADTM, 
PRICE, REVCST,ORDCST , SHOCST , RLAM, SIMRUN 


ieee) SECTION FORMATS THE READ AND WRITE STATEMENTS 


380 FORMAT(1X,"RUN WITH TIME BETWEEN REVIEWS SET AT: "I6,//) 

250 FORMAT(1217,/) 

300 FORMAT(2X,'REQN',1x, ‘HITIME' ,1X,'TOTIME',1X,'REVTIM', 
1X,'REVCNT!',1X,'INVLEV',1X,'INVPOS',1X,'ORDCNT', 
1X, 'SUMORD',1X,'ONORD',1X,'SUMLD',1X, 'TOTSO"') 

FORMAT (1X, 'THEORETICAL:',2X, 'EOQQTY',5X, "MAXINV',5X, 

mere@st' 5X. 'ORDCST' ,5X, "HOLCST' ,5X. "SHOCST' 5x, 

Beevesr',5X,'TOTCST',5X,' TVC *,TX,'AVGSOA') 

FORMAT(12X,10F11.2,/) 
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800 
il 
iL 
900 
810 
SUS 
i 
it 
B75 
2939 
iL 


Oo 
or 
= 


Cae) Co"() (ae 


OO Om © 


60 


65 


68 


FORMAT(1X,'SIMULATION : ',2X,'EOQQTY',5X,'MAXINV' ,5X, 
TMeosm asx ORDeST 5X. 'HOLCST! ,5X,'SHOCST’ , 
See Recor ao. Tones’ asx. TVC '.5x.' SMA 1) 

FORMAT (12X,9F11.2,F11.4,/) 

FORMAT(F5.2,15) 

FORMAT(5X%,'REQNYR',5X,' REVYR',5X,; ORDYR',5X,' SOYR ', 
6X,'SET',7X,'LEADTM' ,5X,'PRICE',5X,'REVCST! ,5X, 
'ORDCST',5X,'SHOCST',5X,' RLAM ',5X,'SIMRUN') 

FORMAT(4F11.2,2111,5F11.2,I11,//) 

FORMAT(1X,'MAX INVENTORY: ',F8.2,3X,'AVG REORDER: ', 
F8.2,3X,'SAFETY STOCK: ',F8.2,3X,'EXPECTED B/O: ', 
F8.4,3X,'GENBRATED B/O: ',F8.4,/) 

FORMAT(1X,'CYCLE DEMAND: ',F6.1,2X,'LEADTIME DEMANDS: ', 
F6.1,2X,'THEORETICAL REVIEWS: ',F6.1,/) 


READ IN THE PERCENTAGE CHANGE IN TIME BETWEEN REVIEWS 
seen CHECK DIGIT FOR END OF RUN 


meebD(5,810) CHANGE, CHECK 
MeCCHECK.GT.0) GO TO 10 
Pepe = DSEED + 233.2D0 


SOP 
END 


Eo, OF PROGRAM 
BEGINNING OF SUBROUTINE SECTION 


fees SUBROUTINE COMPUTES A NEW ORDER. THE QUANTITY ORDER- 
Mmeewiiih BE THE DIFFERENCE BETWEEN THE INVENTORY POSITION 
AND THE MAXIMUM INVENTORY LEVEL 


SUBROUTINE ORDER(MAXINV,ONORD,LEADTM, INVLEV,0C,0L,0Q) 
INTEGER 0C(100) ,OL(100) ,0Q(100) 
INTEGER INVLEV,ONORD, LEADTM 
REAL MAXINV 
NO BO ei Se 

TEC OCG) -EOn@ cen ne. 65 

CONTINUE 
GO TO 68 
ed 
OC(K)=1 
0Q(K) =(IFIX)MAXINV-INVLEV-ONORD 
OL(K) =LEADTM 

RETURN 

END 


bal 





O73) 


THIS ROUTINE UPDATES LEADTIMES AND ADDS TO INVENTORY LEVEL 
BemeeNOLDER QUANTITY LF IT 2S TIME FOR THE ORDER TO COME IN. 
imo sUMS THE REQUISITIONS OUTSTANDING, THEIR QUANTITIES 
AND LEADTIMES 


SuerOUlTINE UPDTCONORD,AITIME,SUMORD,INVLEV,SUMLD,OC,OL, 


he 


80 


70 


Ts 


30 


RETURN 


END 


Ge) o> 


ENTRY 
5s 
0 
mo 

ee 2.0 

P00 

ie 


PhP REP 


OG, INVPOS,ORDCNT , STKOA, AVGSOA) 
INTEGER ONORD,HITIME ,SUMORD,INVLEV,SUMLD 
INTEGER INVPOS,N 
REAL STKOA,AVGSOA 
INTEGER 0C(100) ,OL(100) ,0Q(100) 
DO 80 J=50 

BECOC CE) REOl) Onc ONC) HiT iME 
CONTINUE 
N=51 


DOs? 0 =i 50 
BECO De GE Ve nme On@iore me Hd PiMeE jn N=ul 
CONTINUE 
ir COOWN) ahEn®) Corl a74 
AVGSOA = STKOA/FLOAT(ORDCNT) 
INVLEV=INVLEV+00(N) 
OC(N) =0 
Voc) =o0 
OL(N)=0 


ONORD=0 

SUMOTY=0 

SUMLD=0 

DO 90 J=1,50 
ad 
ONORD=ONORD+0Q(K) 
SUMORD=SUMORD+0C(K) 
SUMLD=SUMLD+OL(X) 

CONTINUE 

INVPOS=INVLEV+ONORD 


fees lo WHERE THE INPUTED CHANGE FROM THE OPTIMAL TIME 
PelweeEN REVIEWS IS ENTERED. USE ZERO TO TERMINATE THE RUNS 
mietnE CHECK FIGURE 


BZ 
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APPENDIX C 


POISSON PROBABILITY DISTRIBUTION PROGRAM 


E 
C THIS PROGRAM PRODUCES A POISSON PROBABILITY DISTRIBUTION 
C BASED ON A DEMAND IN A LEADTIME, OR DEMAND IN A LEADTIME 
C PLUS AND ORDER CYCLE. THE CUMULATIVE DISTRIBUTION 
C FUNCTION (CDF) AND THE COMPLEMENTARY CUMULATIVE DISTRI- 
C BUTION (CCDF) ARE ALSO PRODUCED. 
C 
$JOB 
C 
C DECLARE VARIABLES 
C 
REAL LEADTM,PROB,CDF,CCDF,CCCDF 
INTEGER N 
C 
C INPUT LEADTIME AND CHECK FOR END OF FILE 
C 
eer INT,* ' 
PRINT,! ' 
PRINT, 'PROVIDE NEW LEADTIIE:' 
READ, LEADTM 
IF(LEADTM.LT.O.0) STOP 
é 
C CONPUTE THE POISSON PROBABILITIES 
C 
N = 0 
PROB = EXP(-LEADTM) 
CDF = PROB 


wept = 1.0 - PROB 
eeeDr = 1.0 


C 
meer TITLES FOR DISTRIBUTION AND WRITE PROBABILITIES 
e 
oer, ! 
Pent,’ WN PROB CIDE ENO aD. 
1 CEE PCX > _N) CIDE SOG Se. 
PRINT,! ' 


IF (PROB.GE.0.00001) WRITE (6,10) N,PROB,CDF,CCDF,CCCDF 
N=1l 


Seeels DISTRISUTIONS 
maar LOAT(N).LE.LEADIM) THEN DO 
PROB = PROB * M/FLOAT(N) 
eeecor = 8.0 = ChE 
ener = ChE + PROB 
cepE = 1.0) = CDE 
mCP ROB GE .0 0000 a iertEGo le) N.PROB,CDF , 
Sei. cCCDE 


od ©) 
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1 


ELSE DO 


EiDy FE 
Ste. 0 5 


eS Nae 


PROB = PROB * LEADTM/FLOAT(N) 
IF(PROB.LT.0.00001) GO TO 1 
@coDE.- 1.0. "Cnr 

CDF = CDF + PROB 

CCD = 1,0 Sens 
WRITE(6,10)N,PROB,CDF,CCDF,CCCDF 
N=Ne#il 


Cc 
SeeORMAT STATEMENT 
c 
0 


FORMAT(1X,13,5X,F10.4,F15.4,F19.4,F22.4) 


END 
SENTRY 
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